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ABSTRACT 

In this letter we present a study of the color magnitude relation of 468 early-type galaxies in the 
Virgo Cluster with Sloan Digital Sky Survey imaging data. The analysis of our homogeneous, model- 
independent data set reveals that, in all colors (it — <?, g — r, g — i, i — z) similarly, giant and dwarf 
early-type galaxies follow a continuous color magnitude relation (CMR) that is best described by an 
S-shape. The magnitude range and quality of our data allows us to clearly confirm that the CMR in 
Virgo is not linear. Additionally, we analyze the scatter about the CMR and find that it increases in 
the intermediate-luminosity regime. Nevertheless, despite this observational distinction, we conclude 
from the similarly shaped CMR of semi-analytic model predictions that dwarfs and giants could be 
of the same origin. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: dwarf — galaxies: fundamental 
parameters — galaxies: clusters: individual: (Virgo Cluster) 



1. INTRODUCTION 

That early-type galaxies form a well-defined 
Color Magnitu de R el ation (CMR) was first 
recognized by iBauml (11959D. Ex tensiv e stud- 
ies by IVisvanathan fc Sandagel (|1977l ) and 
ISandage fc Visvanathanl (|1978bl lal) compared the CMR 
of nine galaxy clusters and of field galaxies. They found 
them to be astonishingly similar. In combination with 
other studies comparing the C MR for E and SO, within 
clust ers, groups and the field (jFaberl [l~973t iBower et al.1 
1992) a high de gree of universality was shown for the 
relation. 

The CMR is typically explained by an increase of 
mean stellar metallicity with increasing galaxy mass 
as the dominant effect. The common underlying idea 
is that more massive galaxies have deeper potential 
wells, which can retain metal-enriched stellar ejecta more 
effectively and subsequently recycle the enriched gas 
into new stars ( Kodam a fc Arimotolll997l: iFerreras et al.l 
[1991 [Gailazzi et all 120061 : IChang et al.ll2006f l. Another 
factor could be a variable integrated galactic initial 
mass function, with more massive stars in more mas- 
sive galaxies, and thus a more substantial en richment 
(e.g. iMatteucci et al.lll998t iKoppen et aTll2007f ). 

An alternative explanatio n was given with a ch ange 
of the mean age (see e.g. |Poggianti et al. 200 j, and 
some observational support for it in lRakos fc Schombertl 
I2004T) . The question of which explanation to favor is not 
totally settled yet due to the ambiguity introduced by 
the age-metallicity degeneracy. But data from clusters 
at higher redshift show that the CMR is in place ever 
since, favoring t he mass-metallicity relati on to be the 
dominant effect (|Kodama fc Arimotolll997T ). 

From early on, it was discussed whether the univer- 
sality of the CMR also holds over the whole range of 
galaxy masses, i.e. whether dwarf and giant early-type 
galaxies follow the same CMR. Studies of different clus- 
ters show consistency with one common CMR for dwarfs 



and giants, albeit with a significant increase in the scat- 
ter at low luminosit ies (ISecker et al.l 119971 for Co ma, 
Conselice et al.1 120021 for Perseus^ Karick et al. 2003 and 



Mieske et al.ll2007l for Fornax. ISmith Castelli et al.ll2008l 
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for Antli a and[M isgeld et al.ll2008l for Hydra I). More ex- 
plicitely, I Caldwell (|1983l ) stated that there is a common 
linear relation. But his Fig. 3 might hint at a change 
of slope from h i gh to low luminosities, similar to what 
Ide Vaucouleursl (| 19611 ) suggested. Interestingly, visual 
examination of the diagrams presented by most of the 
above-mentioned studies indicates consistency also with 
a change of slope - yet linear relations we re fitted in most 
cases (see, however. iFerrarese et al.|[2006t see Sect. 5). In 
this letter we revisit the question of the universality of 
the CMR for dwarfs and giants. 

2. SAMPLE SELECTION AND DATA 

The sample is based on the Virgo Cluster Catalog 
(VCC; iBinggeli etafl 119850 . All early-type galaxies 
therein with a certain cluster member status and tub < 
18.0 mag are taken into account, which is the same mag- 
nitude limit up to which the VCC was found to be com- 
plete. This translates into Mb < —13.09 mag with our 
adopt ed distance mo dulus of m-M=31.09 mag (d=16.5 
Mpc, |Mei et a"Lll2007l). We treat uncertain classifications 
as in Uanz fc Liskerl i|2008h and exclude possible irregu- 
lars. Seven further galaxies are exluded due to their low 
S/N in u and z. 

The Sloan Digital Sky Survey (SPSS) Da ta Release 
Five (DR5) (jAdelman-McCarthv et~all [20071 covers all 
but six early- type dwarf galaxies of the VCC. Since the 
quality of sky level subtraction of the SDSS pipeline is 
ins ufficient, we use sk y-subtracted images as provided 
by iLisker et al.l (|2007l ) . based on a careful subtraction 
method. The images were flux-calibrated a nd corrected 
for galactic extinction (Sch legel et al.|[l998T ). 

For each g alaxy, we deter mined a "Petrosian semi- 
major axis" (|Petrosianlll976l hereafter Petrosian SMA, 
a p ), i.e., we use ellipses instead of circles in the cal- 
culat ion of the Petrosian radius (see, e.g., lLotz et aTl 
2004) . The total flux in the r-band was measured within 
a = 2a p , yielding a value for the half-light semimajor 
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Fig. 1. — Color Magnitudes Relations for the colors u — g, g — r, g — i, and i — z. Open circles show dEs with blue cores and filled 
circles all the other galaxies in the sample. Limits for the x-axes are scaled to the mean color of galaxies fainter than M r > — 15 mag and 
brighter than Mr < —21 mag. The gray line indicates the "running histogram" as found in successive magnitude bins with a width of 1 
mag and steps of 0.25 mag, clipped one time at 3o\ We limit the drawing range for the line to the region with at least three galaxies in a 
bin. The white histograms show the distributions in bins of the same width, normalized to the square root of the number of galaxies in the 
bin, shown for every fourth step. The white errorbars indicate typical photometric errors at the respective brightness. In the panels right 
to the color magnitude diagrams a measure for the intrinsic scatter is given by the ratio of the RMS of the scatter and of the photometric 
error in the same running bins (see text). 

axis, CLhl,r,uncorr- This Petrosian aperture still misses 
some flux, w hich is of particula r relevance for the gi- 
ant galaxies (|Truiillo et al.l [200Th . The brightness and 
the half-light SMA were correcte d for this missing flux 
according to iGraham et al.l (|2005f) . 2 Axial ratio and po- 
sition angle were then determined through an isophotal 
fit at a = 2ahi, r - Colors were measured within the el- 
liptical r-band half-light aperture for each filter. Errors 
were estimated from the S/N and calibration uncertain- 
ties (which we estimate to have a relative effect of 0.1 mag 
in each band, which is smaller t han the absolute v alues 
given by SDSS), as described in lLisker et all (|2008l ). 

Our data constitute a very homogeneous set of mea- 
surements for galaxies in one cluster, obtained and re- 
duced with the same instrument, setup, and software. 

2 The correction is computed with the assumption and fitting 
of Sersic profiles as in Uanz fc Lisken pOOl ). But the effect on the 
CMR is small: for a de Vaucouleurs profile it leads to a difference 
of 0.2 mag in brightness and practically no difference in color, since 
the correction to the radius has a minute effect due to the small 
color gradients. 



3. THE CMR OF EARLY- TYPE GALAXIES 

From the five SDSS filter bands we choose four repre- 
sentative colors: u — g, mostly sensitive to age; g — r with 
the highest S/N; g — i with the longest wavelength base- 
line at good S/N; and i — z, mostly sensitive to metal- 
licity. 

First of all, the impression one can get by examin- 
ing just the black points in Fig. 1 is that there is not 
one common linear relation from the faint to the bright 
galaxies. In all colors the overall shape appears more like 
"S" shaped. The brightest (M r < —21) galaxies have al- 
most constant color, i.e. no correlation between color and 
brightness; the very brightest galaxies show a larger scat- 
ter. These were reported before to be morphologically 
different from the other galaxi es in more detailed studies 
of the inner light pr o files fe.g. [Ravindranath et al.ll2001l; 
Truiillo et all 120041: [Ferrarese et al l l2006t lLa uer et al l 
20071 iKormendy et al.ll2009f ). For the remaining galaxies 
several descriptions seem to be plausible, ranging from 
just an offset between two relations with similar slopes 
up to a curved relation. 
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Fig. 2. — Color Magnitude Diagram, comparison to the semi- 
analytic model. Virgo cluster galaxies are shown with black 
symbols and model gala xies from the Numerical Galaxy Catalog 
(Nagashima et~aTl 2005) as gray dots (we selected galaxies with 
B/T > 0.6 and with higher effective surface brightness than 25.5 
mag arcsec - 2 ). Contours were calculated using model galaxy abun- 
dances in bins of 0.5 mag in brightness and 0.02 mag in color. The 
contour levels, relative to the lowest level, are 2, 3.4, 5, 7.5 10, and 
15. 

With the non-linear shape, it seems not very favorable 
to fit a straight line. This would not describe the data 
well, and there is no theoretical prediction what other 
function is expected. So at first, we want to make the 
overall shape more clearly visible, using continous, over- 
lapping magnitude bins, in which mean and scatter are 
calculated. In Fig. 1 these derived relations are shown 
with grey lines. The first impression is confirmed: one 
common linear relation for dwarfs and giants cannot be 
seen. Moreover, the white histograms showing the galaxy 
distributions in the bins are clearly peaked towards the 
bright and the faint end, while they are rather flat at 
intermediate luminosity. 

If the notion of two separate relations with an overlap is 
correct, this should become evident through an increased 
scatter about the relation around the intersection, rela- 
tive to the photometric error. Therefore we calculate 
the RMS of the scatter around the mean in running bins 
(clipping one time at 3cr) and divide it by the RMS of 
the photometric errors 



RMS-ratio 



rms scat _ VEi(p- ( c » 2 



rmserr 



with color c and mean color (c) , averaging over the galax- 
ies in the respective bin. Here we exclude dEs with 
blue cores, since they are known to have different col- 
ors (jLisker et al.ll2008l ). This RMS-ratio should be one 
if the scatter is only due to the measurement errors and 
larger than one if there is an intrinsic scatter. 3 In Fig. 1 
we show the CMRs along with the RMS-ratio p. 

Indeed, the RMS-ratio is enhanced between 
—20 mag < M r < —16 mag, indicating an intrinsi- 
cally increased scatter, which could in principle be 

3 With the binning this is not completely true anymore, since 
flatter relations with the same amount of scatter will show a larger 
RMS-ratio. But if the relation becomes significantly flatter, either 
two relations intersect with an offset anyway, or there is a transition 
region with a flatter CMR; in both cases it is different from the 
relations faint- and brightward. 
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Fig. 3. — Comparison to the data of Fcrrarcsc ct al. (2006). 
Shown are all galaxies in our sample with gray symbols, again 
the blue-core dEs with open circles and the other galaxies with 
filled circles. Galaxies common to both samples are indicated by 
an additional black dot in the middle of the other symbol. The 
axes are changed in comparison to Fig. 1 and colors are measured 
in smaller apertures to most closely resemble their Fig. 123. The 
black lines indicate their fits to their data, the short dashed and 
long dashed line linear fits to a brighter (M 9 < —18) and fainter 
subsample (M g > —18), respectively. The solid line displays the 
parabola they fitted to the whole sample. There is a systematic 
offset of about 0.1 mag tow ards bluer colors of the color derived in 
our study in comparison to Fcrrarcsc ct al. 

explained by a transition between two separate rela- 
tions. But the change is steady. Thus also alternative 
explanations, such as a more varied star formation 
history with decreasing galaxy mass, seem plausible 
- although the decrease in the scatter for the faintest 
dwarfs would still need to be explained. 

One can argue about the significance of the RMS-ratio 
increase for the brightest galaxies, since it is just a hand- 
ful of them - nevertheless, this larger scatter might be re- 
lated to the absence of a well-defined CMR at the bright- 
est magnitudes. 

4. COMPARISON TO SEMI-ANALYTICAL MODEL 



The Numerical Galaxy Catalog of Nagashim a et al] 
(2005) is based on a hig h resolution TV-body simula- 
tion in a ACDM universe (lYahagil 120051) . The dark-halo 
merger trees of the iV-body simulation are taken as in- 
put for a semi-analytical model (SAM) of galaxy for- 
mation (here a mod i fied v ersion of the Mitaka Model, 
iNaeashima k Yoshiil |2004| ). The SAM models physi- 
cal processes of galaxy formation and evolution such as 
gas radiative cooling, star formation, heating by super- 
nova explosions (supernova feedback), mergers of galax- 
ies, population synthesis, and extinction by internal dust 
and intervening HI clouds. In particular, the model 
takes into account the dynamical response to starburst- 
induced gas removal after gas- rich mergers (also for cases 
intermediate between a purely baryonic cloud and a bary- 
onic cloud fully supported by surrounding dark matter 
as in lYoshii fc Arimotol 1 19871 ). The gravitational well 
is shallower for the dwarf galaxies and thus they suffer 
a more substantial gas loss than giants. This process is 
not only important for the galaxy sizes, but also for their 
metal enrichment histories. 

In the comparison of the CMR of the Virgo galaxies 
with the SAM (Fig. 2), the shapes of the distributions 
are fairly similar, and the model CMR is indeed not well 
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represented by a linear relation. In lJanz fc Liskerl (2008) 
we presented an analogous comparison for the scaling re- 
lation of radius and brightness and found that the model 
qualitatively agrees with distribution of the galaxies in 
that diagram. Both the CMR and the relation between 
size and brightness do not show a linear, nor one com- 
mon behaviour for dwarfs and giants. Nevertheless, in 
the framework of the SAM both can be explained by the 
same physical processes, which govern ACDM structure 
formatio n, and thus both c an be of cosmological origin 
(see also IChilingarianl l2009( ) . 

Beside the good agreement in the overall shape, an 
offset is observed. This offset could partly be due to 
uncertainties of the adopted synthetic stellar popula- 
tion model. 4 Furthermore, the relative number of bright 
galaxies exceeds the observed one and the luminosity 
function is clearly different, which could possibly be ex- 
plained with model input physics. 

5. SUMMARY AND DISCUSSION 

In the literature, there is no consistent view of com- 
monness or distinctness of the color magnitude relation 
of dwarf and giant early-type galaxies. Some studies find 
a linear relation over the whole range of brightness, oth- 
ers find a change of slope. 

We studied the CMR of Virgo cluster early-type galax- 
ies, based on the wealth of SDSS imaging data in multi- 
ple filter bands. Our main result is that the dwarfs and 
giants do not follow one common linear relation. The ap- 
pearance in the different colors, from u — g to i — z, is very 
similar, suggesting that age and metallicity go hand in 
hand: the CMR at shorter wavelengths is more sensitive 
to changes in age, while it is more sensitive to changes in 
metallicity at longer wavelengths. Recent studies indeed 
claim that a combination of the effects of age, metallicity 
and also q-enhancement shape the CMR (|K5ppen et al.l 
[20071: iLisker fc Hadl20f)l . 

The m ost direct comparison with other data is possi- 
ble with iFerrarese et all (|2006T ). who studied early- type 
galaxies in the Virgo Cluster with the Hubble Space Tele- 
scope's Advanced Camera for Surveys (ACS) in the ACS 
Virgo Cluster Survey. In their Fig. 123 they show the 
CMR in g — z, using smaller apertures for the color mea- 
surements. For comparison we repeat the plot and indi- 
cate galaxies common to both samples with black dots 
in the symbols (Fig. 3). We also show their fits (to their 
data): two linear fits to a brighter and fainter subsam- 
ple and one parabola for the whole sample. While a 
qualitatively good agreement is found, the small offset 

4 The offset might here also be influenced to some degree by the 
transformations to Vega magnitudes in the Johnson-Cousins filter 
system. In a diagram of i — z vs. Mi (for which no transformations 
are needed) the observed galaxies lie bluewards of the model ones 
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to bluer colors of our data is likely explained by a com- 
bination of the effect of different apertures, the atmo- 
sphere, and diff erent physical filters. While we agree with 
IFerrarese et al.l that the CMR is not linear over the whole 
range, the faint dwarfs in our sample, which they do not 
reach (cf. Fig. 3), do not lie on the extrapolation of their 
parabolic fit, but rather define an overall "S"-like shape. 
Furthermore, some of the galaxies that might drag their 
CMR to the blue side, at their faint end, display blue 
cores from recent star formation activity (|Lisker et al.l 
l2006h . which we excluded. 

We found an increase of the scatter about the CMR 
at intermediate brightness, visible in Fig. 1 through a 
broadening of the color distributions within the respec- 
tive magnitude bins. This confirms the result that dwarfs 
and giants do not share one common linear CMR. For 
the brighter and fainter galaxies the scatter is closer to 
t he photometric er rors. 

iBernardi et~aT1 (|2005h concluded that the CMR is 
a result of two more fundamental relations: the 
Faber-Jackson relation and a relation between color 
and velocity dispersion. Given the slope change 
of the Faber-Jackson relati o n from giants t o dw arfs 
(|Matkovic fc Guzman! 120051 : Ide Riicke et~aT1 l2005f) a 
change of slope of the CMR would actually be expected. 

It is important to emphasize that the curved shape 
of the CMR, and the possible existence of two adja- 
cent relations does not necessarily imply different for- 
mation scenarios of dwarfs and giants. This is in ac- 
cordance with previous claims of no distinction be- 
tween them (|Graham fc Guzman! 120031: iGavazzi et all 
120051: IFerrarese et al.ll2006l:lMisgeld et al.ll2008f ). Instead, 
from the good qualitative agreement of the observed 
CMR shape with a semi-analytic model, we conclude 
that a common origin within the framework of ACDM 
structure formation appears plausible for both giant and 
dwarf early types. A more detailed multiparameter com- 
parison of dwarfs and giants will be communicated in a 
forthcoming paper. 
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with an smaller absolute offset than in B — V vs. Mb, but the 
scatter of the observed galaxies is larger due to the lower S/N in 
the z-band. 
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